arXiv:1502.03864v2 [physics.atom-ph] 29Jun2015 


Stability of a trapped atom clock on a chip 


R. Szmuk, 1 V. Dugrain, 2 W. Maineult, 1 J. Reicliel, 2 and P. Rosenbusch 1 

1 LNE-SYRTE, Observatoire de Paris, PSL Research University, 

CNRS, Sorbonne Universites, UPMC Univ. Paris 06, 

61 avenue de I’Observatoire, 75014 Paris , France 
2 Laboratoire Kastler Brossel, ENS, UPMC, CNRS, 24 rue Lhomond, 75005 Paris, France 

(Dated: June 30, 2015) 

We present a compact atomic clock interrogating ultracold 8 ' Rb magnetically trapped on an atom 
chip. Very long coherence times sustained by spin self-rephasing allow us to interrogate the atomic 
transition with 85% contrast at 5 s Ramsey time. The clock exhibits a fractional frequency stability 
of 5.8 x 10 -13 at 1 s and is likely to integrate into the 10 -15 range in less than a day. A detailed 
analysis of 7 noise sources explains the measured frequency stability. Fluctuations in the atom 
temperature (0.4 nK shot-to-shot) and in the offset magnetic field (5 x 1CU 6 relative fluctuations 
shot-to-shot) are the main noise sources together with the local oscillator, which is degraded by the 
30% duty cycle. The analysis suggests technical improvements to be implemented in a future second 
generation set-up. The results demonstrate the remarkable degree of technical control that can be 
reached in an atom chip experiment. 


I. INTRODUCTION 

Atomic clocks are behind many everyday tasks and 
numerous fundamental science tests. Their performance 
has made a big leap through the discovery of laser cool¬ 
ing m giving the ability to control the atom position 
on the mm scale. It has led to the development of atomic 
fountain clocks ms which have reached a stability lim¬ 
ited only by fundamental physics properties, i.e. quan¬ 
tum projection noise and Fourier-limited linewidth [6|. 
While these laboratory-size set-ups are today’s primary 
standards, mobile applications such as telecommunica¬ 
tion, satellite-aided navigation [7| or spacecraft naviga¬ 
tion [5 call for smaller instruments with litre-scale vol¬ 
ume. In this context, it is natural to consider trapped 
atoms. The trap overcomes gravity and thermal expan¬ 
sion and thereby enables further gain on the interroga¬ 
tion time. It makes interrogation time independent of 
apparatus size. Typical storage times of neutral atoms 
range from a few seconds to minutes 0 OH. Thus a 
trapped atom clock with long interrogation times could 
measure energy differences in the mHz range in one single 
shot. Hence, if trap-induced fluctuations can be kept low, 
trapped atoms could not only define time with this reso¬ 
lution, but could also be adapted to measure other physi¬ 
cal quantities like electromagnetic fields, accelerations or 
rotations with very high sensitivity. A founding step to¬ 
wards very long interrogation of trapped neutral atoms 
was made in our group through the discovery of spin 
self-rephasing El which sustains several tens of seconds 
coherence time 5M3|. This rivals trapped ion clocks, 
the best of which has shown 65 s interrogation time and 
a stability of 2 10 -14 at 1 s [THUS]- It is to be compared 
to compact clocks using thermal vapour and buffer gas 
ushu] or laser cooled atoms [T5H22I- Among these the 
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record stability is 1.6 x 10 13 at 1 s m- Clocks with 
neutral atoms trapped in an optical lattice have reached 
impressive stabilities down to the 10~ 18 range but 

their interrogation time is so far limited by the local os¬ 
cillator. Research into making such clocks transportable 
is on-going 05J155]. We describe the realization of a com¬ 
pact clock using neutral atoms trapped on an atom chip 
and analyze trap-induced fluctuations. 

Our ’’trapped atom clock on a chip” (TACC) employs 
laser cooling and evaporative cooling to reach ultra-cold 
temperatures where neutral atoms can be held in a mag¬ 
netic trap. Realising a 5 s Ramsey time, we obtain 
100 mHz linewidth and 85% contrast on the hyperfine 
transition of 87 Rb. We measure the fractional frequency 
stability as 5.8 x 10 -13 t - 1 / 2 . It is reproduced by analyz¬ 
ing several noise contributions, in particular atom num¬ 
ber, temperature and magnetic field fluctuations. The 
compact set-up is realized through the atom chip tech¬ 
nology mi which builds on the vast knowledge of micro¬ 
fabrication. The use of atom chips is also widespread for 
the study of Bose-Einstein condensates 0123, degener¬ 
ate Fermi gases |29| and gases in low dimensions ;3Ql [31] ■ 
Other experiments strive for the realization of quantum 
information processors The high sensitivity and 

micron-scale position control have been used for prob¬ 
ing static magnetic [35] and electric [36] fields as well 
as microwaves m ■ Creating atom interferometers pIBj 
on atom chips is equally appealing. Here, an on-chip 
high stability atomic clock not only provides an excellent 
candidate for mobile timing applications, it also takes a 
pioneering role among this broad range of atom chip ex¬ 
periments, demonstrating that experimental parameters 
can be mastered to the fundamental physics limit. 

This paper is organised as follows: we first describe the 
atomic levels and the experimental set-up. Then we give 
the evaluation of the clock stability and an analysis of all 
major noise sources. 
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II. ATOMIC LEVELS 

We interrogate the hyperfine transition of s 'Rb (fig¬ 
ure [I]). A two photon drive couples the magnetically 
trappable states | 1 ) = |F=l,mj? = — 1 ) and | 2 ) 

I = 2 ,rriF = 1 ), whose transition frequency exhibits a 
minimum at a magnetic field near B m « 3.229 G p9ll40j . 
This 2 nd order dependence strongly reduces the clock fre¬ 
quency sensitivity to magnetic field fluctuations. It as¬ 
sures that atoms with different trajectories within the 
trap still experience similar Zeeman shifts. Furthermore, 
by tuning the offset magnetic held, the inhomogeneity 
from the negative collisional shift [Hi] can be compen¬ 
sated to give a quasi position-invariant overall shift Bl¬ 
under these conditions of strongly reduced inhomogene¬ 
ity we have shown that spin self-rephasing can overcome 
dephasing and that coherence times of 58 ± 12 s m can 
be reached. It confirms the possibility to create a high 
stability clock [42]. 

~ 2 MHz 

F = 2 


F = 1 


m F =— 2-1 0 1 2 

FIG. 1: Level scheme of the s 'Rb ground state. The 
two magnetically trappable clock states 
\F = l,rriF = —1) and |F = 2 ,itif = 1) are coupled via 
a 2 -photon, microwave and radiofrequency transition, 
where the microwave is tuned 500 kHz above the 
| 1 , — 1 ) —> | 2 , 0 ) transition. 



III. EXPERIMENTAL SET-UP 

The experimental set-up, details of which are given 
in [43] . is similar to compact atom chip experiments re¬ 
ported previously |44}|45]. All experimental steps, laser 
cooling, evaporative cooling, interrogation and detection 
take place in a ~ (5 cm ) 3 glass cell where one cell wall 
is replaced by the atom chip (figure [5]). In this first- 
generation set-up, a 25 1/s ion pump is connected via 
standard vacuum components. It evacuates the cell to 
a pressure of ~ 1 x 10 ” 9 mbar. The cell is surrounded 
by a 10 x 10 x 15 cm 3 cage of Helmholtz coils. A 30 cm 
diameter optical table holds the coil cage as well as all 
beam expanders necessary for cooling and detection and 
is surrounded by two layers of magnetic shielding. 


The timing sequence (table [I]) starts with a mirror 
MOT [44] loading ~ 3 x 10 6 atoms in ~ 7 s from the 
background vapor. The MOT magnetic field is gener¬ 
ated by one of the coils and a U-shaped copper struc¬ 
ture placed behind the atom chip [H^. Compressing 
the MOT followed by 5 ms optical molasses cools the 
atoms to ~ 20 /iK. The cloud is then optically pumped 
to the |1) state and transferred to the magnetic trap. It is 
gradually compressed to perform RF evaporation, which 
takes ~ 3 s. A 0.7 s decompression ramp transfers the 
atoms to the final interrogation trap with trap frequen¬ 
cies (u> x ,u> y ,ui z ) =27r x (2.7,92,74) Hz located 350 /mr 
below the surface. It is formed by two currents on the 
chip and two currents in two pairs of Helmholtz coils. The 
currents are supplied by homebuilt current supplies with 
relative stability < 10 ~ 5 at 3 A 03- The final atom num¬ 
ber is N = 2 — 4 x 10 4 and their temperature T ~ 80 nK. 
The density is thus with n ~ 1.5 x 10 11 atoms/cm 3 so 
low that the onset of Bose-Einstein condensation would 
occur at 5 nK. With kpT/Hui x , yiZ > 20 the ensem¬ 
ble can be treated by the Maxwell Boltzmann distri¬ 
bution. The trap lifetime 7 -1 = 6.9 s is limited by 
background gas collisions. The clock transition is in¬ 
terrogated via two-photon (microwave + radiofrequency) 
coupling, where the microwave is detuned 500 kHz above 
the |1) to |F = 2 ,mp = 0) transition (figure]!]). The mi¬ 
crowave is coupled to a three-wire coplanar waveguide 
on the atom chip f!3] US] . The interaction of the atoms 
with the waveguide evanescent field allows to reach sin¬ 
gle photon Rabi frequencies of a few kHz with moderate 
power ~ 0 dBm. Since the microwave is not radiated, 
interference from reflections, that can lead to field-zeros 
and time varying phase at the atom position, is avoided. 
Thereby, the waveguide avoids the use of a bulky mi¬ 
crowave cavity. The microwave signal of fixed frequency 
vmw ~ 6.8 GHz is generated by a homebuilt synthesiser 
)48l which multiplies a 100 MHz reference signal derived 
from an active hydrogen maser [49] to the microwave fre¬ 
quency without degradation of the maser phase noise. 
The actual phase noise is detailed in section VB The 
RF signal of variable frequency vpp ~ 2 MHz comes 
from a commercial DDS which supplies a ’’standard” 
wire parallel to the waveguide. The two-photon Rabi 
frequency is about f2 = 3.2 Hz making a 7 t/2 pulse last 
t p = 77.65 ms 27 t/uj z . The pulse duration is chosen so 
that any Rabi frequency inhomogeneity, which was char¬ 
acterised in EDI, is averaged out and Rabi oscillations 
show 99.5% contrast. Two pulses enclose a Ramsey time 
of Tp = 5 s. Detection is performed via absorption imag¬ 
ing. A strongly saturating beam crosses the atom cloud 
and is imaged onto a back illuminated, high quantum ef¬ 
ficiency CCD camera with frame transfer (Andor iKon 
M 934-BRDD). 20 /rs illumination without and with re¬ 
pump light, 5.5 ms and 8.5 ms after trap release, probes 
the F = 2 and F = 1 atoms independently. Between 
these two, a transverse laser beam blows away the F = 2 
atoms. Numerical frame re-composition generates the re¬ 
spective reference images and largely reduces the effect 
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of optical fringes m- Calculation of the optical density 
and correction for the high saturation [52] give access to 
the atom column density. The so found 2D atom dis¬ 
tributions are fitted by Gaussians to extract the number 
of atoms in each state 1 Vi, 2- The transition probability 
is calculated as P = 7V 2 /(1Vi + 7V 2 ) accounting for total 
atom number fluctuations. The actual detection noise is 


discussed in section |V A| The total time of one experi¬ 
mental cycle is T c = 16 s. 

a b 



FIG. 2: (Color online) (a) The implemented atom chip. 
One identifies the Z-shaped coplanar waveguide which 
serves for atom trapping and transport of the 
microwave interrogation signal. The outer dimensions 
are 38 x 45.5 mm 2 , (b) The chip constitutes one of the 
facets of the vacuum cell facilitating electrical contact. 
The cell is surrounded by a 10 x 10 x 15 cm 3 cage of 
Helmholtz coils and a 30 cm diameter optical table 
holding all optical beam expanders. The cell is 
evacuated via standard UHV equipment. 


Operation 

Duration 

MOT 

6.85 s 

compressed MOT 

20 ms 

optical molasses 

5 ms 

optical pumping 

1 ms 

magnetic trapping and compression 

230 ms 

RF evaporation 

3 s 

magnetic decompression 

700 ms 

first Ramsey pulse 

77.65 ms 

Ramsey time 

5 s 

second Ramsey pulse 

77.65 ms 

time of flight (1) , 2)) 

(5.5 , 8.5) ms 

detection 

20 y s 


TABLE I: Timing sequence of one experimental cycle. 
The total cycle time is 16 s. 


IV. STABILITY MEASUREMENT 

Prior to any stability measurement we record the typ¬ 
ical Ramsey fringes. We repeat the experimental cycle 
while scanning v^o = vmw + v rf over ~ 3 fringes. Do¬ 
ing so for various Ramsey times Tr allows to identify the 
central fringe corresponding to the atomic frequency is at . 



FIG. 3: Typical Ramsey fringes recorded at Tr = 5 s 
while scanning the local oscillator detuning. Each point 
corresponds to a single experimental realization. One 
identifies the Fourier limited linewdith of 100 mHz and 
the very good contrast of 85%. 


Figure [3] shows typical fringes for Tr = 5 s, where each 
point is a single shot. One recognises the Fourier limited 
linewidth of 100 mHz equivalent to ~ 10 11 quality factor. 
The 85% contrast is remarkable. A sinusoidal fit gives the 
slope at the fringe half-height dP/dv = 13.4/Hz, which is 
used in the following stability evaluation to convert the 
detected transition probability into frequency. 

Evaluation of the clock stability implies repeating the 
experimental cycle several thousand times. The clock is 
free-running, i.e. we measure the transition probability 
at each cycle, but we do not feedback to the interroga¬ 
tion frequency vlo- Only an alternation in successive 
shots from a small fixed negative to positive detuning, 
A mo ^/(27 t) = ±50 mHz, probes the left and right half¬ 
height of the central fringe. The difference in P between 
two shots gives the variation of the central frequency in¬ 
dependent from long-term detection or microwave power 
drifts. In the longest run, we have repeated the frequency 
measurement over 18 hours. 

The measured frequency data is traced in figure [4] ver¬ 
sus time. Besides shot-to-shot fluctuations one identifies 
significant long-term variations. Correction of the data 
with the atom number, by a procedure we will detail in 
section V C 1 results in substantial improvement. We 
analyze the data by the Allan standard deviation which 
is defined as |M1 


IL/2'J-I 

a y( T ) =2 (^+! _ dkf ■ (1) 

k=l 

Here L is the total number of data points and the yk are 
averages over packets of 2 l successive data points with l S 
{0,1, ...|k>g 2 LJ} and r = 2 l T c . Figure [ 5 ] shows the Allan 
standard deviation of the uncorrected and corrected data. 
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For 0 < l < 9 the points and their errorbars are plotted 
as calculated with the software ”Stable32” [54]. This 
software uses equation [l] to find the points. The error 
bars are calculated as the 5% - 95% confidence interval 
based on the appropriate x' 2 distribution. The software 
stops output at l = Llog 2 L\ — 2 since there are too few 
differences yk+i ~Vk to give a statistical errorbar. Instead 
we directly plot all differences for l = 10 and 11. 

The Allan standard deviation shows the significant 
improvement brought by the atom number correction. 
The uncorrected data starts at r = T c = 16 s with 
a y = 1.9 x 1CT 13 shot-to-shot. For the TV-corrected data, 
the shot-to-shot stability is a y = 1.5 x 10~ 13 . Up to 
r ~ 100 s the corrected frequency fluctuations follow a 
white noise behaviour of ctj,(t) = 5.8 x 10 -13 t -1 / 2 . At 
r ss 1000 s, the fluctuations are above the r _1//2 be¬ 
haviour but decrease again at r > 5000 s. For r > 10 4 s, 
3 of the 4 individual differences are below 10” 14 . This 
lets us expect that a longer stability evaluation would in¬ 
deed confirm a stability in the 10 -15 range with sufficient 
statistical significance. The ’’shoulder” above the white 
noise behaviour is characteristic for an oscillation at a 
few 10 3 s half-period. Indeed, this oscillation is visible in 
the raw data in figure [4] Its cause is yet to be identified 
through simultaneous tracking of many experimental pa¬ 
rameters - a task which goes beyond the scope of this 
paper. 

Table [TT] gives a list of identified shot-to-shot fluc¬ 
tuations that contribute to the clock frequency noise. 
Treating them as statistically independent and summing 
their squares gives a fractional frequency fluctuation of 
1.5 x 10 -13 shot-to-shot or 6.0 x 10 -13 at 1 s, correspond¬ 
ing to the measured stability. We have thus identified 
all major noise sources building a solid basis for future 
improvements. In the following we discuss each noise 
contribution in detail. 


Relative frequency stability (10 16 ) 

shot-to-shot 

at 1 s 

measured, without correction 

2.0 

7.2 

measured, after TV correction 

1.5 

5.8 

atom temperature 

1.0 

3.9 

magnetic field 

0.7 

2.6 

local oscillator 

0.7 

2.7 

quantum projection 

0.4 

1.5 

TV correction 

0.4 

1.5 

atom loss 

0.3 

1.1 

detection 

0.3 

1.1 

total estimate 

1.5 

6.0 


TABLE II: List of identified contributions to the clock 
(in)stability. Atom temperature fluctuations dominate 
followed by magnetic field fluctuations and local oscilla¬ 
tor noise. The quadratic sum of all contributions explains 
the measured stability. 


• Clock data 



§. With atom number correction 
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FIG. 4: (Color online) Relative frequency deviation 
when repeating the clock measurement over 18 h, (top) 
raw data, (bottom) after correction by the 
simultaneously detected total atom number. The blue 
dots represent single shots, red dots show an average of 
10 shots. 



FIG. 5: (Color online) Allan standard deviation of the 
measured clock frequency with (blue circles) and 
without (red diamonds) atom number correction. For 
integration times smaller than 10 4 s, the points and 
errorbars are calculated using the software ” Stable32”. 

Above 10 4 s, the individual differences between 
successive packets of 1024 and 2048 measurements are 
given. The TV-corrected data follows initially 
5.8 x 10 -13 r -1 / 2 (blue dashed line). The quantum 
projection noise and the local oscillator noise are given 
for reference. 


V. NOISE ANALYSIS 


In a passive atomic clock, an electromagnetic signal 
generated by an external local oscillator (LO) interacts 
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with an atomic transition. The atomic transition fre¬ 
quency v a t is probed by means of spectroscopy. The de¬ 
tected atomic excitation probability P is either used to 
correct the LO on-line such that vlo = Fat, or, as applied 
here, the LO is left free-running and the measured dif¬ 
ferences [vlo ~ y at){t) are recorded for post-treatment. 
The so calibrated LO signal is the useful clock output. 

When concerned with the stability of the output fre¬ 
quency, we have to analyze the noise of each element 
within this feed-back loop, i.e. 

A. noise from imperfect detection, 

B. folded-in fluctuations of the LO frequency known as 
Dick effect, 

C. fluctuations of the atomic transition frequency in¬ 
duced by interactions with the environment or between 
the atoms. 

We begin by describing the most intuitive contribution 
(A. detection noise) and finish by the most subtle (C. 
fluctuations of the atomic frequency). 


A. Detection and quantum projection noise 

The clock frequency is deduced from absorption imag¬ 
ing the atoms in each clock state as described in section 
|III| N-\ and N 2 are obtained by fitting Gaussians to the 
atom distribution, considering a square region-of-interest 
of ~ 3 x 3 cloud widths. 

Photon shot noise and optical fringes may lead to atom 
number fluctuations of standard deviation o^et- These 
fluctuations add to the true atom number. Analyzing 
blank images, we confirm that a^ et increases as the num¬ 
ber of pixels in the region-of-interest and that optical 
fringes have been efficiently suppressed m- This scal¬ 
ing has led to the choice of short times-of-flight where 
the atoms occupy fewer pixels [523- Supposing the same 
crdet for both states, we find for the transition probability 
noise a Pjdet = a det /(V2N) with N = N x + N 2 . 

Another P degradation may occur if the Rabi fre¬ 
quency of the first pulse fluctuates or if the detection 
efficiency varies between the |1) and |2) detection. The 
latter may arise from fluctuations of the detection laser 
frequency on the timescale of the 3 ms difference in tirne- 
of-flight. Both fluctuations induce a direct error ap P f + if 
on P independent from the atom number. 

Quantum projection noise is a third cause for fluctua¬ 
tions in P. This fundamental noise arises from the fact 
that the detection projects the atomic superposition state 
onto the base states. Before detection, the atom is in a 
near-equal superposition of |1) and |2). The projection 
then can result in either base state with equal proba¬ 
bility giving ctqpn = 1/2 for one atom. Running the 
clock with N (non-entangled) atoms is equivalent to N 
successive measurement resulting in up^qpn = 1 /(2y/~N) 
shot-to-shot. 

We quantify the above three noise types from an inde¬ 
pendent measurement: Only the first ir/2 pulse is applied 


and P is immediately detected. The measurement is re¬ 
peated for various atom numbers and crp(N) is extracted. 
Figure [6] shows the measured a P shot-to-shot versus N. 
Considering the noise sources as statistically indepen¬ 
dent, we fit the data by <jp = cr^ et /(2A 2 ) + 1/(4AT) + 
a p RF+if an( i Vdet = 59 atoms and (J P ,rf+i f < 10” 4 . 
ddet is equivalent to an average of 2.2 atoms/pixel for 
our very typical absorption imaging system. The low 
ciRF+if proves an excellent passive microwave power sta¬ 
bility < 2.5 x 10 -4 , which may be of use in other exper¬ 
iments, in particular microwave dressing ESI EH- 

During the stability measurement of figure [4] about 
20 000 atoms are detected, which is equivalent to 
&y,QPN = 0.4 x 10 -13 shot-to-shot. The detection region- 
of-interest is slightly bigger than for the above characteri¬ 
sation, so that adet = 69 atoms, corresponding to a y det = 
a d et{PatN)- l \dP/dv r 1 = ^ = 0.3 x 10- 13 

shot-to-shot. In both we have used dP/dv as measured 
in figure [3] 



FIG. 6: Characterization of the detection noise. Only a 
single 7 t/ 2 pulse is applied and P = N 2 /N detected. 

The shot-to-shot Allan deviation is plotted as a 
function of the total atom number. We fit the data with 
the quadratic sum of the detection noise adet/{V^N), 
the quantum projection noise 1/(2 y/N) and the Rabi 
frequency and laser frequency noise ap,Rf+if. The fit 
gives <j de t = 59 atoms and crp^pf+if < 10 -4 . 


B. Local oscillator noise 

The experimental cycle probes v at — vpo only during 
the Ramsey time. Atom preparation and detection cause 
dead time. Repeating the experimental cycle then consti¬ 
tutes periodic sampling of the LO frequency and its fluc¬ 
tuations. This, as well-known from numerical data acqui¬ 
sition, leads to aliasing. It folds high Fourier frequency 
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LO noise close to multiples of the sampling frequency 
1 /T c back to low frequency variations, which degrade the 
clock stability. Thus even high Fourier frequency noise 
can degrade the clock signal. The degradation is all the 
more important as the dead time is long and the duty cy¬ 
cle d = Tji/T c is low. This stability degradation a Vt Dick 
is known as the Dick effect [53]. It is best calculated us¬ 
ing the sensitivity function g(t) [53]: during dead-time, 
g = 0 whereas during Tr, when the atomic coherence 
\tp) = (11) + |2))/v^2 is fully established g = 1. Dur¬ 

ing the first Ramsey pulse, when the coherence builds 
up, g increases as sin f 'It for a square pulse and decreases 
symmetrically for the second pulse ED]. Then the inter¬ 
rogation outcome is 


I-TJ2 K t(t) - v L o(t)) g(t) dt 
f-TJ 2 9(t) ^ 


( 2 ) 


with 


g(t) 


{ a sin D(Tr /2 + t p + t) 
a sin Qt p 

asinf2(^p + t p — t) 

0 


Is < + < Is 4 . r 

2 — 1 — 2 r /p 

otherwise. 


( 3 ) 

Typically flr p = 7 r /2 and, for operation at the fringe half¬ 
height, a = sinA mo( jTfj = 1. Because of the periodicity 
of repeated clock measurements, it is convenient to work 
in Fourier space with 


gi = ^r[ g(t)cos(2nlt/T c )dt. (4) 

J-c J -T c /2 


Using the power spectral density of the LO frequency 
noise S^(f), the contribution to the clock stability be¬ 
comes the quadratic sum over all harmonics m 

I °o , \ 2 

? £ [fj S>(l/Ta). (5) 


Here we have assumed v at constant in time; its fluctu¬ 
ations are treated in the next section. The coefficients 
idi/g o ) 2 are shown as points in figure [ 8 ] for our condi¬ 
tions. The weight of the first few harmonics is clearly the 
strongest, rapidly decaying over 6 decades in the range 
1/T C « 0.1 Hz to 1/Vp ss 10 Hz. Above ~ 10 Hz the gi 
are negligible. 

To measure S^{f) we divide our LO into two princi¬ 
pal components: the 100 MHz reference signal derived 
from the hydrogen maser and the frequency multipli¬ 
cation chain generating the 6.8 GHz interrogation sig¬ 
nal. We characterise each independently by measuring 
the phase noise spectrum S^(f). The fractional fre¬ 
quency noise S£ (/) is obtain from simple differentiation 
as S£(f) = f 2 S</,(/) / v 2 mw [59]. The frequency noise of 
the RF signal can be neglected as its relative contribution 
is 3 orders of magnitude smaller. 


We characterize the frequency multiplication chain 
by comparing it to a second similar model also con¬ 
structed in-house. The two chains are locked to a com¬ 
mon 100 MHz reference and their phase difference at 
6.834 GHz is measured as DC signal using a phase detec¬ 
tor (Miteq DB0218LW2) and a FFT spectrum analyzer 
(SRS760). The measured S${f) is divided by 2 assuming 
equal noise contributions from the two chains. It is shown 
in figure 7} It features a 1// behaviour up to / = 10 Hz 
and reaches a phase flicker floor of —115 dB rad 2 /Hz at 
1 kHz. The peak at / = 200 Hz is due to the phase lock 
of a 100 MHz quartz inside the chain to the reference 
signal. As we will see in the following, its contribution 
to the Dick effect is negligible. 

The 100 MHz reference signal is generated by a 
100 MHz quartz locked to a 5 MHz quartz locked with 
40 mHz bandwidth to an active hydrogen maser (VCH- 
1003M). We measure this reference signal against a 
100 MHz signal derived from a cryogenic sapphire oscilla¬ 
tor (CSO) [511 162] . Now the mixer is M/A-COM PD-121. 
The CSO is itself locked to the reference signal but with 
a time constant of ~ 1000 s [B3j. This being much longer 
than our cycle time, we can, for our purposes, consider 
the two as free running. The CSO is known from prior 
analysis [M| to be at least 10 dB lower in phase noise than 
the reference signal for Fourier frequencies higher than 
0.1 Hz. Thus the measured noise can be attributed to the 
reference signal for the region of the spectrum / > 1/T C 
where our clock is sensitive. The phase noise spectrum 
is shown in figure [7J For comparison it was scaled to 6.8 
GHz by adding 37 dB. Several maxima characteristic of 
the several phase locks in the systems can be identified. 
At low Fourier frequencies, the reference signal noise is 
clearly above the chain noise. For all frequencies, both 
are well above the noise floor of our measurement sys¬ 
tem. The noise of the reference signal being dominant in 
the range 1/T C to l/r p , where our clock is sensitive, we 
neglect the chain noise in the following. 

Using equation [5] we estimate the Dick effect contri¬ 
bution as CTy t Dick = 2.7 x 10” 13 t - 1 / 2 . This represents 
the second biggest contribution to the noise budget (ta¬ 
ble [TT|. It is due to the important dead time and the long 
cycle time which folds-in the LO noise spectrum where it 
is strongest. Improvement is possible, first of all, through 
reduction of the dead time which is currently dominated 
by the ~ 7 s MOT loading phase and the 3 s evaporative 
cooling. Options for faster loading include pre-cooling 
in a 2D MOT [65] or a single-cell fast pressure modula¬ 
tion [BB]. Utilization of a better local oscillator like the 
cryogenic sapphire oscillator seems obvious but defies the 
compact design. Alternatively, generation of low phase 
noise microwaves from an ultra-stable laser and femtosec¬ 
ond comb has been demonstrated by several groups iZb 
1551 and on-going projects aim at miniaturisation of such 
systems [70]. If a quartz local oscillator remains the pre¬ 
ferred choice, possibly motivated by cost, one long Ram¬ 
sey time must be divided into several short interrogation 
intervals interlaced by non-destructive detection D3HZ3]. 
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FIG. 7: (Color online) Phase noise power spectral 
density of the local oscillator. The frequency 
multiplication chain and the 100 MHz reference signal 
are characterised separately. The beat between two 
quasi identical chains is performed at 6.8 GHz (red). 
The beat of the reference signal against a cryogenic 
sapphire oscillator is taken at 100 MHz and scaled to 
6.8 GHz (black). The noise of the reference signal 
dominates in the low frequency part, where our clock is 
sensitive. Both results are above the intrinsic noise of 
the measurement system (blue). 
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FIG. 8: (Color online) (black) Same data as figure [7] 
now expressed as fractional frequency fluctuations 
S y = (red) Fourier coefficients of the 

sensitivity function ( gi/go ) 2 for our conditions 
(T c = 16 s, Tr = 5 s and t p = 77.65 ms). Multiplication 
of the two gives the stability degradation known as Dick 

effect. 


FIG. 9: Correlation between the detected atom number 
and the clock frequency for the data of figure [4] Fitting 
with a linear regression gives k = —2.70(7) /iHz/atom, 
which allows to correct the clock frequency at each shot 
and yields substantial stability improvement. 


quency. We begin by atom number fluctuations. Due to 
the trap confinement and the ultra-cold temperature, the 
atom density is 4 orders of magnitude higher than what 
is typically found in a fountain clock. Thus the effect of 
atom-atom interactions on the atomic frequency must be 
taken into account even though 87 Rb presents a substan¬ 
tially lower collisional shift than the standard 133 Cs. In¬ 
deed, when plotting the measured clock frequency against 
the detected atom number N = N\ + N%, which fluctu¬ 
ates by 2-3% shot-to-shot, we find a strong correlation 
(figure [9]) . The distribution is compatible with a linear 
fit with slope k = —2.70(7) /.iHz/atom. In order to con¬ 
firm this value with a theoretical estimate we use the 
mean field approach and the s-wave scattering lengths 
a,ij which depend on the atomic states only 001 


A vc{r) = — n(r) [(a 2 2 - an) + (2a i2 - an - 022 )$] ■ 

( 6 ) 


n(r) is the position dependent density and an = 
100.44ao, a 2 2 = 95.47ao, a -12 = 98.09ao are the scattering 
lengths with a 0 = 0.529 x 10 -10 m ^Dj. We assume per¬ 
fect 7 t/ 2 pulses and so 9 = (Ni — = 0. Integrating 

over the Maxwell-Boltzmann density distribution we get 


C. Fluctuations of the atomic frequency 

1. Atom number fluctuation 


A vc = N x 


-fr{au - a, 22 )\/mu} x u}yuj z 
4(7t fc^T) 3 / 2 


(7) 


Having characterised the fluctuations of the LO fre¬ 
quency, we now turn to fluctuations of the atomic fre- 


We must consider that the atom number decays during 
the Tr = 5 s since the trap life time is 7 -1 = 6.9 s. We 




























replace N by its temporal average 

rT R 


N = — I 

Tr Jo 


= N ., 
= N f 


Nie-^dt 
1 - e~ lTR 


lT R 

e l T R _ 1 


( 8 ) 


7 t r 
1.47 N f 


where Ni and Nf are the initial and final atom numbers, 
respectively. Note that Nf is the detected atom number. 
Using T = 80 nK, which is compatible with an inde¬ 
pendent measurement, we recover the experimental colli- 
sional shift of k = —2.7 /iHz/(detected atom). It is equiv¬ 
alent to an overall collisional shift of A vc = —54 mHz 
for Nf = 20 000. 

Using k and the number of atoms detected at each 
shot we can correct the clock frequency for fluctuations. 
The corrected frequency is given in figure [4] showing a 
noticeable improvement in the short-term and long-term 
stability. The Allan deviation indicates a clock stabil¬ 
ity of 5.8 x 10 -13 t -1 / 2 at short term as compared to 
7.2 x 10~ 13 t -1 / 2 for the uncorrected data. At long term 
the improvement is even more pronounced. This demon¬ 
strates the efficiency of the IV-correction. Furthermore, 
the experimentally found k shows perfect agreement with 
our theoretical prediction so that the theoretical coeffi¬ 
cient can in future be used from the first shot on without 
the need for post-treatment. 

While we have demonstrated the efficiency of the atom 
number correction, the procedure has imperfections for 
two reasons: The first, of technical origin, are fluctu¬ 
ations in the atom number detectivity as evaluated in 
section IV A1 The second arises from the fact that atom 
loss from the trap is a statistical process. For the first, we 

get <Jy 

, correction — V%\k\crdet/vat = 0.4 x 10 13 shot-to- 
shot. This value is well below the measured clock stabil¬ 
ity, but may become important when other noise sources 
are eliminated. It can be improved by reducing the atom 
density and thus k or by better detection, in particular at 
shorter time-of-ffight where the camera region-of-interest 
can be smaller. The second cause, the statistical nature 
of atom loss, translates into fluctuations that in principle 
cannot be corrected. The final atom number Nf at the 
end of the Ramsey time is known from the detection, but 
the initial atom number TV,; can only be retraced with a 
statistical error. To estimate this contribution we first 
consider the decay from the initial atom number N. At 
time t, the probability for a given atom to still be trapped 
is e _7t and the probability to have left the trap is 1—e _7t . 
Given N, the probability p to have N t atoms at t is pro¬ 
portional to e~ Ntryt (1 — e d) Nz Nt and to the number 
of possible combinations: 


p{N t given Ni) 


m 


s—Ntjt 


N t \(Ni — N t )\ 


(l-e- 7t ) 
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The sum of this binomial distribution over all 0 < N t < 
N is by definition normalised. We are interested in the 
opposite case: since we detect Nf at t = T R , we search 
the probability of N t given Nf. 


p{N t given Nf) = 


AN t \ 

N f \(N t — N f )\ 


X e ~ N fl( T R- t ) (1 _ e -'Y(TR-t)\N t -N f 

( 10 ) 


The combinatorics are as in equation [9] when replacing 
Nt —► Nf and N —> Nt, but now normalisation sums 
over 0 < Nt < oo. Here it is convenient to approximate 
the binomial distribution by the normal distribution 

p{Nt given Nf) « -£= e ~ (N f - N t e~^R~^ /(2r,) (n) 

v2tt?7 

with r) = Nte~' 1 ^ TR ~ t \\ — e~ 1< ' TR ~ t ' > ) and hence A = 
e -~t(T R -t)' Then, the mean of N t is 

(N t ) = (N f + l)e 7(Tj? -* ) - 1 w N f e' l(TR - t) (12) 

and its statistical error 


<j Nt = ^(1 - e 7 ( T «-*))(2- (N f +2)eA TR -t)) 
« sjN f (eA TR -t) - l ) e ATR-t)_ 


(13) 


Setting t = 0, we get <t/v. = 210. Integrating ajy t over 
Tr gives <t]v = 113 « a Nil 2 and a frequency fluctuation 
of (J y ,Losses = 0.3 x 10~ 13 shot-to-shot. This can be im¬ 
proved by increasing the trap lifetime well beyond the 
Ramsey time, which for our set-up implies better vac¬ 
uum with lower background pressure. Alternatively one 
can perform a non-destructive measurement of the ini¬ 
tial atom number m- Assuming an error of 80 atoms 
on such a detection would decrease the frequency noise 
to (7 y .losses = 0.1 x 10” 13 shot-to-shot. 


2. Magnetic field and atom temperature fluctuations 

We have analyzed the effect of atom number fluctua¬ 
tions. Two other parameters strongly affect the atomic 
frequency: the atom temperature and the magnetic field. 
We show that their influence can be evaluated by mea¬ 
suring the clock stability for different magnetic fields at 
the trap center. We begin by modelling the dependence 
of the clock frequency. 

Our clock operates near the magic field B m « 3.229 G 
for which the transition frequency has a minimum of - 
4497.31 Hz with respect to the field free transition, 

A Mn = b(B(r)-B m ) 2 (14) 

with b « 431 Hz/G 2 . For atoms trapped in a harmonic 
potential in the presence of gravity, the Zeeman shift be¬ 
comes position dependent 

A v B (r) = (ulx 2 + uly 2 + lv 2 z z 2 - 2gz + 6B-) 

fi B V y m / 

(15) 


( 9 ) 
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with SB = B{f = 0) — B m and g the gravitational accel¬ 
eration [41]. Using the Maxwell-Boltzmann distribution 
the ensemble averaged Zeeman shift is 


Az's = -2~( 
Ms 


b Ag 2 mk B T 


+ 15 klT 2 


(16) 


6g B SBk B T + 5B 2 g %). 


Differentiation with respect to SB leads to the effective 
magic field 


SB o B 


-3 k B T 
Mb 


(17) 


where the ensemble averaged frequency is independent 
from magnetic field fluctuations. For T = 80 nK, 
SBq = — 3.6 mG whose absolute value almost coincides 
with the magnetic field inhomogeneity across the cloud 
/-—\ 1/2 

f B{f) 2 -B 2 ) = V6k B T/g B = 2.92 mG. SB^ is 

close to the field of maximum contrast SBq ss —40 mG 
such that the fringe contrast is still 85% (figure [l0|) . 

If SB ^ SBq is chosen the clock frequency fluctuations 
due to magnetic field fluctuations are 


a y,B 


2b 


V a t 


SBq 


SB\ <j b . 


(18) 


We will use this dependence to measure a B . 

Temperature fluctuations affect the range of magnetic 
fields probed by the atoms and the atom density, i.e. the 
collisional shift. Differentiation of both with respect to 
temperature also leads to an extremum, where the clock 
frequency is insensitive to temperature fluctuations. The 
extremum puts a concurrent condition on the magnetic 
field with 


5B% = 


15 k B T+^ 

3 g B 

h{an - a 2 2 )(e 7Tfl - 1 )y/mN f g B ui x u)yU z 


16n 3 / 2 b(k B T) 5 / 2r yT B 


(19) 


For our conditions, SBq = —3.6 mG and SBq = —79 mG 
are not identical but close and centered around SBq . We 
will see in the following that a compromise can be found 
where the combined effect of magnetic field and temper¬ 
ature fluctuations is minimised. A ’’doubly magic” field 
can not be found as always SBq < SBq , but lower T re¬ 
duces their difference. If SB ^ SBq is chosen, the clock 
frequency fluctuations due to temperature fluctuations 
are 


a = ^b_\ 5b t - 5 B\a T ( 20 ) 

Mb v a t. 

thus varying SB allows to measure <jt , too. 

We determine a B and experimentally by repeat¬ 
ing several stability measurements for different SB over 



FIG. 10: Fringe contrast (top) and differential Zeeman 
shift (bottom) of the clock frequency with respect to 
the frequency minimum for various magnetic fields. 
SB = 0 indicates the magic field of 3.229 G. The 
contrast maximum is offset by —40 mG. 


a range of 200 mG where the contrast is above 70%. The 
shot-to-shot stability is shown in figure |TT] One identi¬ 
fies a clear minimum of the instability at SB ss —35 mG, 
which coincides with SBq and is a compromise between 
the two optimal points SBq and SBq. This means, 
that both magnetic field and temperature fluctuations 
are present with roughly equal weight. We model the 
data with a quadratic sum of all so far discussed noise 
sources. Most of them give a constant offset; the slight 
variation due to the contrast variation shown in figure [10| 
is negligible. a ViB and a Vt r are fitted by adjusting a B 
and We find shot-to-shot temperature fluctuations 
of ctt = 0.44 nK or 0.55% relative to 80 nK. The shot- 
to-shot magnetic field fluctuations are a B = 16 g G or 
5 x 10~ 6 in relative units. The values demonstrate our 
exceptional control of the experimental apparatus. Be¬ 
cause the ambient magnetic field varies by < 10 mG and 
the lowest magnetic shielding factor is 3950, we attribute 
a B to the instability of our current supplies. Indeed, it is 
compatible with the measured relative current stability 
m- The atom temperature fluctuations are small com¬ 
pared to a typical experiment using evaporative cooling. 
This may again be due to the exceptional magnetic field 
stability, since the atom temperature is determined by 
the magnetic field at the trap bottom during evapora¬ 
tion and the subsequent opening of the magnetic trap. 
At all stages, the current control is the most crucial. 
Using equations [18] and [20j the temperature and mag¬ 
netic field fluctuations translate into a frequency noise 
of er^T = 1-0 x 1CU 13 and a y , B = 0.7 x 10 -13 shot-to- 
shot, respectively. The comparison in table [TT] shows, 
that these are the main sources of frequency instability 
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together with the Dick effect. Therefore, improving the 
magnetic field and temperature noise is of paramount 
importance. The atom temperature can in principle be 
extracted from the absorption images, which we take at 
each shot. Analysis of the data set of figure [4] gives 
shot-to-shot fluctuations of <jt/T = 2 — 4%, which is 
much bigger than the 0.55% deduced above. We there¬ 
fore conclude that the determination of the cloud width is 
overshadowed by a significant statistical error. Neverthe¬ 
less, it needs to be investigated, whether better detection 
and/or imaging at long time-of-flight, may reduce this 
error. The magnetic field stability may be improved by 
refined power supplies, the use of multi-wire traps 
microwave dressing m or ultimately the use of atom 
chips with permanent magnetic material )76H78| . If the 
magnetic field fluctuations can be reduced, the tempera¬ 
ture fluctuations may also reduce. Small ctb would also 
allow to operate nearer to SBq . 



FIG. 11: (Color online) Shot-to-shot clock stability for 
various magnetic fields. Error bars are smaller than the 
point size. One observes a clear optimum at 
SB = —35 mG. Fitting with the quadratic sum of all 
identified noise contributions allows to quantify the 
atom temperature fluctuations (0.4 nK shot-to-shot) 
and magnetic field fluctuations (16 nG shot-to-shot). 
The individual contributions are shown as dashed lines. 

Two sweet spots exist where the temperature 
dependence and the magnetic field dependence vanish. 


VI. CONCLUSION 


We have built and characterised a compact atomic 
clock using magnetically trapped atoms on an atom chip. 
The clock stability reaches 5.8 x 10 -13 at 1 s and is 
likely to integrate into the 10 -15 range in less than a day. 
This is similar to the performance of the best compact 
atomic microwave clocks under development. It further¬ 
more demonstrates the high degree of technical control 
that can be reached with atom chip experiments. After 
correction for atom number fluctuations, variations of the 
atom temperature and magnetic field are the dominant 
causes of the clock instability together with the local os¬ 
cillator noise. The magnetic field stability may be im¬ 
proved by additional current sensing and feedback and 
ultimately by the use of permanent magnetic materials. 
This would allow to operate nearer to the second sweet 
spot where the clock frequency is independent from tem¬ 
perature fluctuations. The local oscillator noise takes an 
important role, because the clock duty cycle is < 30%. 
We are now in the process of designing a second version 
of this clock, incorporating fast atom loading and non¬ 
destructive atom detection. We thereby expect to reduce 
several noise contributions to below 1 x 10~ 13 t -1 / 1 2 3 4 5 . 
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